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Lie group methods are used for the study of various issues related to symmetries and 
exact solutions of the barotropic vorticity equation. The Lie symmetries of the barotropic 
vorticity equations on the /- and /?-planes, as well as on the sphere in rotating and rest 
reference frames, are determined. A symmetry background for reducing the rotating 
reference frame to the rest frame is presented. The one- and two-dimensional inequivalent 
subalgcbras of the Lie invariance algebras of both equations are exhaustively classified and 
then used to compute invariant solutions of the vorticity equations. This provides large 
classes of exact solutions, which include both Rossby and Rossby-Haurwitz waves as special 
cases. We also discuss the possibility of partial invariance for the /3-plane equation, thereby 
further extending the family of its exact solutions. This is done in a more systematic and 
complete way than previously available in literature. 



1 Introduction 

The governing equations of geophysical fluid dynamics are mainly nonlinear partial differential 
equations (PDE). Since there is no general theory available for solving such equations, it is 
known to be very difficult to systematically construct their exact solutions. In meteorology, this 
problem is usually overcome by solving the governing equations numerically. However, as models 
become more sophisticated, it may be difficult to directly evaluate the quality of these numerical 
results. Moreover, it is dissatisfactory to rely solely on numerical modeling when studying the 
physics of the atmosphere. It is thus to be expected that exact solutions can both enhance our 
understanding of atmospheric processes and provide consistency tests for numerical models. 

The classical method of reduction of PDEs by using its Lie symmetries |19t [20] and the exten- 
sion to partially invariant solutions [20] provides a manageable way to systematically construct 
exact solutions. It is the goal of this paper to carry out a comprehensive symmetry investigation 
of the barotropic vorticity equation both on the /3-plane and on the sphere. Although there 
are already a number of works on the /3-plane equation [U [H [I3l [Ml [HI [16], none of them 
gives a systematic and complete symmetry analysis. In Ref. [3] the classification of inequiva- 
lent subalgebras is done only for the one-dimensional case. The symmetry properties were used 
in Ref. [5] in order to obtain new solutions from the known ones. In a recent paper [13], the 
procedure of group-invariant reduction is done without reference to the algebraic aspects of the 
classification problem. Consequently, some of the reductions presented in Ref. [13] are overly 
complicated, and hence in some cases, these authors were only able to obtain some particular 
solutions (most notably of the well-known Rossby wave class). This reveals that the vorticity 
equation has classes of completely integrable reduced PDEs, as shown in the present paper. 
Finally, Refs. \15\ I16| (see also Ref. \14:\ pp. 221-225]) also contain only a nonsystematic list of 
some group-invariant solutions. To the best of our knowledge, the spherical equation has not 
been investigated in light of its symmetries at all so far. 



We divide this paper into two main parts: the first deahng with the symmetry analysis 
of the equation on the /3-plane and the second considering the spherical version. For both 
equations, we determine the maximal Lie invariance algebras and classify their one- and two- 
dimensional subalgebras. Based on this classification, we give a complete list of group-invariant 
reduced equations and then demonstrate that Rossby (Rossby-Haurwitz) waves can be realized 
as group- invariant solutions of the barotropic vorticity equation on the plane (the sphere). Also, 
by means of algebraic inspection of the Lie symmetry algebras, it is shown that for the spherical 
equation there is no need to consider rotation of the Earth. Finally, some examples for partially 
invariant solutions will be given for the /3-plane equation. 

2 The /3-plane equation 

This section contains the classical symmetry analysis of the barotropic vorticity equation on the 
/3-plane (/3BVE). 

2.1 The model 

Assuming the two-dimensional velocity field v to be nondivergent, it is possible to cast the Euler 
equations of an ideal fluid in a rotating reference frame as the conservation law of absolute 
vorticity rj = C + /> where C = k • (V x v) is the vertical component of the vorticity vector 
(relative vorticity) and / denotes the vertical Coriolis parameter, which depends only on y. In 
what follows, we approximate / by its truncated Taylor series, / = /o + d//dy|oy =: fo + Pu, 
which leads to the /3-plane approximation [12]. The Euler equations can then be equivalently 
written as the /?BVE 

Ci + i'xCy - IpyCx + P^Px = 0, (1) 

where we have used the shorthand notation Q = d(/dt, etc. The stream function = ip[t,x,y) 
generates a nondivergent flow. It is related to the vorticity by means of the Laplacian, 

Rescaling allows us to set /3 = 1, but for physical reasons this is not desired here. 

2.2 The symmetries 

The barotropic vorticity equation can be considered as a submodel of the ideal Euler equations, 
which have been thoroughly investigated in light of their symmetries (see, e.g., Refs. [Tl I23j). 

Nevertheless, it is instructive to consider the symmetries of the barotropic vorticity equation 
separately to work out the peculiarities of large scale, two-dimensional fluid dynamics. It is 
quite common for different models of incompressible fluids that they admit infinite dimensional 
maximal Lie invariance algebras of a special structure. Techniques for handling such infinite- 
dimensional Lie algebras in order to solve hydrodynamic equations are given, e.g., in Ref. [?]• 
To the best of our knowledge, the symmetry algebra of the barotropic vorticity equation in the 
regular case /3 7^ 0, as well as some exact solutions, was first computed in Refs. [15l |T6] (see 
also Refs. [lllM]). The fact that the singular case /? = admits nontrivial symmetries has been 
known for a long time [3]. The corresponding maximal Lie symmetry algebra was rigorously 
calculated in Ref. [2] (see also Ref. [T]). It is significantly larger than for the regular case of /?. 

We recomputed the symmetry algebras for our purposes and checked them with the computer 
algebra programs MuLie [TO] and DESOLV [6]. In the singular case /? = corresponding to 
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dynamics on the /-plane, the vorticity equation admits the infinite dimensional Lie symmetry 
algebra Bq^ with the basis generators 

Vi = tdt - i)d^ 
J = -ydx + xdy 
dt 

X{f) = f{t)d, - f'{t)yd^, 

where h, f and g run through the set of real-valued time-dependent functions. The shorthand 
notation for partial derivatives, e.g., dt = d/dt is used. The physical significance of these gener- 
ators is as follows: Pi and are scaling operators, dt generates time translations, and and 
J'* correspond to rotations and time-dependent rotations in the horizontal plane. The operators 
y{h) and X{f) are the infinitesimals of generalized transformations on a time-dependent mov- 
ing coordinate system in the y- and x-directions, respectively. The generator 2{g) represents 
gauging the stream function. 

Likewise, in the case /? 7^ 0, eqn. ([T]) admits an infinite-dimensional Lie symmetry algebra 
Bf^, which is a subalgebra of Bq^. The basis generators of Bf^ are 

V = Vi-V2 = tdt-xd,-ydy-3tl;d^, dt, y{l) = dy, X{f), Z{g). 

The physical importance of this algebra is obvious from the /? = case. As shown in the next 
section, this remarkable difference between the cases of vanishing and nonvanishing /3 has no 
counterpart in spherical coordinates. 

For the sake of completeness, we mention that the vorticity equation admits the discrete sym- 
metries (t, X, y, ip) I— > [—t, —X, y, ip) and {t, x, y, ■0) 1— > {t, x, — y, —^p) as well as their composition 
and their compositions with continuous symmetries. 



T>2 = 2^1)8^ + xdx + ydy 

J* = -tyd^ + txdy + ^{x'^ + y'^)d^ 

^{g) = g{t)di, 

y{h) = h{t)dy + h'{t)xd^ 



2.3 Classification of subalgebras 

For an efficient and systematic computation of invariant solutions of PDEs, it is crucial to classify 
their Lie symmetry subalgebras. This is done upon using the adjoint action of a Lie group on its 
Lie algebra, which allows to determine the simplest representatives of equivalent subalgebras. 
The adjoint action of exp(ev) on wq is defined as the Lie series, 

°° e'^ 

w(e) = Ad(exp(ev))wo := V -r{v",wo}, 

n=0 

where we introduced a shorthand notation for nested commutators: {v", wq} := wq, {v", wq} := 
(— l)"[v, {v"^-'^, Wo}]. Alternatively, the adjoint representation can also be calculated by inte- 
grating the initial value problem 

dvif(e) r / N n / \ 

— = w(e), V , w(0) = Wq. 

de 

The nonidentical adjoint actions with basis elements of the algebra Bf^ are exhausted by the 
following list: 

Ad(e^^*)P = V-edt Ad{e'^^''^)V = V + eX{v + tv') 

M{e^^y )V = V + edy Ad(e^^(") )V = V + eZ{'iu + tu') 

kd{e''^)dt = e'dt Ad(e"^("))5t = dt + e^(n') 

Ad(e"'^(^))5t = dt + eX{v') Ad{e'^)dy = e-'dy 
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Ad(e='^(^))5j, = dy- eZ{v') kd{e''^)X{v) = X{v), v = e-^v{e-H) 

Ad{e'^^)X{v) = X{v), v = v{t- e) Ad(e^^-)^(/) = X{f) + Z{f') 

Ad(e"^)Z('u) = Z(n), u = e-^'u{e-H) M{e^^')Z{u) = Z{u), u = u{t-e). 

They are used subsequently to classify the one- and two-dimensional subalgebras of Bf^. 

The approach to the classification of one- dimensional subalgebras is fairly inductive: one takes 
the most general form of an infinitesimal generator of B^, v = a£)T> + afdt + aydy-\-X{f)-\-Z{g), 
and subsequently tries to simplify it using adjoint actions and scaling by a nonvanishing constant 
multiplier [TO]. This is done under additional assumptions on the constants a£), at, and ay and 
the functions f{t) and g{t). Finally, the optimal set of conjugacy inequivalent one-dimensional 
subalgebras of B^ reads 

{V), {dt + cdy), {dy + X{f)), {X{f)+Z{g)), (2) 

where c G {0, ±1}. By means of using the discrete symmetry (t, x, y, i— > (t, —y, —ijj) we can 
further assume c G {0, 1}. Moreover, due to adjoint actions, there are additional equivalences 
inside the third and fourth cases. In the third case, we can apply the adjoint actions Ad{e^''^) to 
rescale the argument t and the function / and Ad(e^^') to shift the argument t of the function /, 
respectively. In the fourth class, the additional equivalences are generated by Ad(e^^'), Ad(e^'^), 
Ad(e^^^) and scaling the basis elements. So, the subalgebras {X{f) + Z{g)) and {X[f) + Z{g)) 
are equivalent if and only if f{t) = af{e^'^t + ei), g{t) = ag{ef'^t + ei) + e^f^ef^t + ei) for some 
constants ei, £2) ^s, and a, where a 7^ 0. 

The procedure for the classification of two-dimensional subalgebras is quite the same as in the 
one-dimensional case: one takes the most general form of two (linearly independent) infinitesimal 
generators = d'-j^V -\- a\dt + aydy + X{f^) + Z{g^) with i = 1,2 and tries to simultaneously cast 
them in a simpler form. This can be done by taking their nondegenerate linear combinations 
and acting on them by adjoint actions under different assumptions on the constants o^, a], and 
Oy and/or functions /* and g*. Since the generators form a subalgebra, one additionally has 
to ensure that the commutator of and lies in their span. This places further restrictions 
on both the constants a^, aj, and a'y and on the functions 

Since the classification of two-dimensional subalgebras is somewhat lengthy, we only give the 
final result here. The list of inequivalent algebras reads 

{V,dt), {V, dy + aX{l)), {V, A:{\tn + cZ{\tr')), {V,Z{\tr^)), 
{dt + bdy, Xie"') + Z{{abt + c)e"*)> , {dt + bdy, Z{{abt + c)e"*)> , 
{dy + X{f^), X{1) + Z{g^)) , {dy + X{f^), Z{g^)) , {X{f^) + Z{g^), X{f) + Z{g')) , 

where o, 6,and c are arbitrary constants and /* = and g^ = g^{t), i = 1,2, are arbitrary 
functions of time. In the last Lie subalgebra, the pairs of functions {f^,g^) and have to 

be linearly independent. Again, within of the above classes there is an additional equivalence 
due to adjoint actions and changes of the basis. 

Computing the associated group-invariant solutions based on the above classification deter- 
mines those cases that have to be considered for receiving the complete set of inequivalent 
solutions. 

2.4 Group-invariant reduction with one-dimensional subalgebras 

In what follows, we give the complete list of reduced equations obtained by imposing invariance 
under the one-parameter groups associated with the Lie algebras ([2]). In all four cases, p,q 
denote the new independent variables, while v = v(p,q) is the new dependent variable. In the 
last case, an ansatz for ip exists only for the values of t, where / 7^ 0. 
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{dt + cdy) = v 

— nin -1— n 7/1 — 1) 11) -\— fii) — n 


11) ' — 


p = X, q = y — ct 

11 —1— 1) 

PP ^ 


3 


{dy + x{f)) ^ = ^_ 1/y 

(1 + p)wq + 2//'u; + f3vp - f" = 


= 0, 


P = x- fy, q = t 
w := Vpp 


4 


{X{f)+Z{g)) ^ = y-txy+jx 
u;g+(f-fp)u;p + /3(f-fp) 


= 0, 


p = y, q = t 

w := Vpp 



We now discuss some properties and/or explicit solutions derived from the different cases 
considered above. 

Case 1. The reduced equation admits the maximal Lie invariance algebra {dp + qd^, d^). The 
basis operators dp + qd^ and 9„ are induced by the operators X{t~^) and Z{t~^) from B^, 
respectively. Hence, there are no hidden symmetries related to this reduction. 

Case 2. The reduced equation admits the maximal Lie invariance algebra {dp, dq, dv) or {pdp + 
qdq + 3vdv, dp, dq, d^) if c = 1 or c = 0, respectively. Any operator from this algebra is obviously 
induced by an operator from B'^. For the basis operators, we have the following correspondence: 

dp <— A'(l), dq <— dy, dv <— Z{1), pdp + qdq + Svd^ < V. Again, we have no related hidden 

symmetries. 

As a result, further Lie reductions in both of the above cases give no new solutions in com- 
parison to Lie reduction with respect to two-dimensional subalgebras. 

Case 3 admits an exhaustive description for its general solution. Integrating once with respect 
to p and using a change of coordinates yields 

Vpq + Pv = 0, (3) 

where 

/" , h , {{1 + P)f"y . f dq 

" = rT7^-7^ + ^ + — ' '' = jTTp^ ^ = ^' 

and h = h{q) is an arbitrary smooth function ol q = t. Eqn. ([3]) is the one-dimensional Klein- 
Gordon equation (presented in the light-cone variables). For a list of exact solutions of this 
equation see, e.g. |22| . It is straightforward to recover the famous Rossby wave solution upon 
using a harmonic ansatz for v and choosing either / = (one-dimensional Rossby waves) or 
/ = const (two-dimensional Rossby waves) [3]. 

Case 4 is completely integrable by quadratures. First, we determine w by solving the char- 
acteristic system. Afterwards, we integrate twice with respect to p to determine v. Finally, 
substituting the expression so obtained for v into the ansatz for ijj, we arrive at the correspond- 
ing group-invariant solution 

= ^F{e) - \[iy^ + h'y + h'- ^xy + fx, 

where /, g, h} , and hP are arbitrary smooth functions of t, F is an arbitrary smooth function 
oi 9 = fy — J gdt. The functions and can be set equal to by symmetry transformations 
generated by an operator of the form X{f) + Z{g). 
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2.5 Group-invariant reduction with two-dimensional subalgebras 



Having considered reduction with one-dimensional subalgebras, it is not overly difficult to in- 
vestigate reduction with two-dimensional subalgebras as well. Namely, the general solutions 
of cases 3 and 4 from section 12.41 are completely described. That is, it is not necessary to 
consider reduction with two-dimensional subalgebras containing the generators dy + X{f) and 
X{f) + Z[g). Moreover, since all algebras containing Z{g) cannot be used for a classical Lie 
reduction, the number of cases that need to be examined reduces to: 

1 i'Djdt) ^ = -sj (2;^ + y^Y" v{(f) 99 = arctan I 

v{w + PsiuLf)^ - lv^{w + l3sin(p) =0, w := v^^ + 9v 

The first of the above equations implies the following functional relation between w and v: 
w + Psimp = cqv^. If Co = 0, the second equation can be easily integrated with respect to v. 
This leads to the invariant solution 

= ci(x^ - 3y^)x + C2(3x2 - y'^)y - ^(x^ + y^)y 

of ([1]). In the case cq 7^ 0, we find particular solutions of the second equation which give rise to 
the invariant solutions 

/ /^^ 2 , 2n| • 3 A 4^ y\ I P I 2 , 2-.^ ■ 3/^1 . y , ^ 

w = — (x + w ) 2 sm - arctan — , w = (x +y )^ sm - arctan — it — 

2.6 Partially invariant solutions 

For a system with at least two dependent variables, it is possible to determine partially invariant 
solutions j20]. The construction of partially invariant solutions has already been extensively 
considered in hydrodynamics p!} [HI II H 117^ 123). In this part, we compute some partially invariant 
solutions for the /3BVE. First of all, it is noted that any single equation can be split into a system 
of multiple equations in various ways introducing a new dependent variable for each additional 
equation desired. We consider the /?BVE as the system of two PDEs 

Ct + i'xCy - i'yCx + Pi^x = 0, ^ = + ^Pyy, (4) 

where both and C are treated as dependent variables. The splitting of eqn. ([T]) into system 
is quite natural since both tp and C have an obvious physical importance. Of course, it is not 
unique. Another natural splitting is given by the system in terms of the usual velocity variables 
together with the condition of vanishing divergence. However, here we will not pursue any other 
splittings further. 

It is an important property of the chosen splitting that the maximal Lie invariance algebra 
of ([4]) is isomorphic to the algebra B"^. More precisely, every operator from B"^ is a prolongation 
of an operator from B"^. This is why for the construction of partially invariant solutions we can 
use the lists of subalgebras obtained above. 

As an example for a partially invariant solution, we use the subalgebra {X{l),Z{g)). Due to 
the generator Z{g), we cannot make an ansatz for However, we can make an ansatz for C and 
because of the generator dx, we have ip = ip{t,x,y) and ( = C{t,y). Therefore, ^ is reduced to 
Ct + '4'xiCy + /3) = 0, C = V'xx + fpyy Introducing the absolute vorticity f] = + j3y and setting 
■0 = X, y) + C{t, y) with C,yy = C, we find rjt + ^x^ly = 0, and '^xx + '^yy = 0. If rjy = 0, we 
have rjt = and, consequently, rj = const. The stream function constructed in this way then 
reads as ^ 

^ = ^ _ i3yS + rjy'^ 
2 
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where '^{t,x,y) is an arbitrary solution of the Laplace equation + "^yy = 0- 
In case rjy ^ 0, we find that the stream function has the form 

= j^Hu^) - l^y' - '-^^x + + f, 

where uj = g^y + and g^,g^ , f^, and are functions of t. 

To present one more example of a partially invariant solution, we take the subalgebra 
{dy,Z{g)). Similar to the previous case, we now have ( = Cit^x) and ^p = 'ip{t,x,y). Then, 
^ is reduced to Q — ipyCx + (^'ipx = and C = V'xx + tpyy Introducing ip = (p{t, x, y) + cj(t, x), 
(^xx = Cj this set of equations yields 

(t>xx + (Pyy = 0, axxt " (^xxx(t>y + P{(f>x + CTx) = 0. (5) 

We now have to distinguish different cases for the integration of this system. 
!• CTxxx = 0. In this instance, the solution for the stream function reads as 

where x'^ are smooth functions of t. The functions a-nd x^ can be set equal to 

by symmetry transformations generated by an operator of the form ^(/) + Z[g). 

2. (Jxxx 7^ 0. We set (f) = H — (3~^axt — f and substitute into the second equation of ([5]), which 
yields a characteristic system for H. Solving this system, we find that H = H(t,rj), where 
V = CTxx + f3y is again the absolute vorticity. From the first equation of ([5]) , we then derive 

H,,i{axxxf + (3^) + H,axxxy - ^ - a.. = 0. (6) 

If we fix X in the above equation, we can write h?'{t)Hj^^ + h}-{t)H^ + h?{t) = 0. We now have to 
distinguish whether there are two independent equations of this type or only one. 

In case of two equations we have i^^^ = and consequently H = a{t)r] + 7(t). Substituting 
this into ([6]), solving the resulting PDE and transforming back to the original variables, we find 

= x) + a{t)(3y, aTi^x — P~^'Sxt-^ P~^S{t) — S = 0. By means of symmetry transformations 
generated by ^(/) and 2(g) we can set a = 6 = and again arrive at the Klein-Gordon 
equation. This illustrates the fact that in some cases the ansatz for a partially invariant solution 
effectively reduces to a usual group-invariant reduction. 

If we only have one independent equation in H, h? ^ 0, and the equation 

H,,{{axxxf + + H,axxxx - ^ - a.. = Kh'm,, + h\t)H, + h\t)l 

where A = \{t,x), holds identically in H. Splitting this equation with respect to H leads to the 
three equations 

(Cxxx)^ + 0^ = A/l^, (Jxxxx = ^h^, —(3 ^O'xxxt — (^xx = A/i*^. 

Since ^ 0, we can express A from the first equation. Provided that = 0, we integrate the 
second equation to find a = a^{t)x^ + a'^{t)x'^ + a^{t)x + cr^{t). Inserting this expression in the 
third equation and splitting with respect to x then yields = 0, i.e., axxx = 0, contradicting 
the initial assumption for this case. For ^ 0, we integrate the second equation once with 
respect to x and then substitute the resulting expression for axx into the third equation. This 
leads to a contradiction in the system constructed by splitting with respect to x, and hence no 
solution is obtained also under the assumption ^ 0. 
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3 The spherical equation 



3.1 The model 

The barotropic vorticity equation on the sphere (sBVE) is given by (e.g. Ref. |21j ^ 



G + (V'aCm 



0, C: = 



1 

i?2 



1 



1 



rV'AA+((l-M'W^ 



(7) 



where ip is the (spherical) stream function and C the (spherical) vorticity. They are related 
through the Laplacian on the sphere. Instead of using the latitude <p as an independent vari- 
able, in practice, it is convenient to rather use /i = siny?. The value of /x ranges from —1 (South 
Pole) to 1 (North Pole). By A we denote the longitude, R is the mean radius of the Earth and 
Q the absolute value of the Earth's angular rotation vector. 



3.2 The symmetries 

We aim to start with ([7]) in a nonrotating reference frame ($7 = 0). Note in passing that it is 
possible to scale the radius R of the Earth to 1 by including R in the stream function via setting 

i> = i)/R^. 

The corresponding Lie symmetry algebra is infinite dimensional and a suitable basis is 
provided by 



V -- 

J2 



tdt - ipd^, 
sin A 



dt, Z{g) = g{t)d^, 



dx + - A*^ cos A5^ , Ja 



Ji = dx, 
cos A 



-.dx - sin A5^. 



(8) 



As for the physical meaning of these basis elements, we find that T) is the generator of scaling 
in t and and dt corresponds to time translations. The generators Ji, i = 1,2,3, correspond 
to rotations in angular coordinates. This follows since they satisfy the commutation relations 
of the Lie algebra so(3), [J7i,v7j] = Ylk=i^ijk'^k, where i,j = 1,2,3 and Eijk is the Levi-Civita 
symbol. 2(g) again represents gauging of the stream function. 

The algebra Sq° has the structure of so(3) ® (02 € {2{g))), where 02 = {T^,dt) is the two- 
dimensional nonabelian algebra and {Z{g)) is an infinite-dimensional Abelian ideal in 5q°. 

Now turning to the rotating case (0 7^ 0). Eqn. ([Tj) admits the infinite-dimensional Lie 
invariance algebra 



V = tdt — [ip — ^fJ')d^ — ^Itdx, 



dt, 



Z{g) =g{t)d^,, 



Ji = dx, 



J2 = ^ — )_ ^ dx H ((1 - ^ + ^d^) , 



cos(A + ^ 
^ „ ^A 



sin(A + ^t) 



11^ 



((1 - ^2)5^ + ^Q^ 



The physical interpretation of the basis elements is obvious from those of the case i7 = 0. 
Moreover, straightforward calculation shows that both Lie symmetry algebras and S"^ are 
isomorphic and can be mapped to each other by means of the change in the coordinates. 



?=t, /i = /i, A = A + Ot, ip = tp — Qfj,. 



(9) 



Furthermore, it is possible to transform ([7]) into the corresponding equation in the rest frame 
($7 = 0) upon using ([9|). This recovers, in a systematic way, the transformation used by Platzman 
|21j to reduce the spherical vorticity equation to a reference frame with zero angular momentum. 
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Note that this mapping is possible due to the special form of the Laplacian in spherical 
coordinates. In particular, it is impossible to obtain a similar result for the vorticity equation in 
Cartesian coordinates since in this case the respective Lie symmetry algebras are nonisomorphic. 
Consequently, no transformation can be found that maps the vorticity equation on the /3-plane to 
the vorticity equation on the /-plane. This indicates that the traditional /3-plane approximation 
significantly distorts the geometry of the more natural spherical vorticity dynamics. 

Again there are two discrete symmetries, given by {t, A, fi, V) ^ —A, fi, ^) and {t, A, fj,, V') ^ 
{t, A, —jj,, —ijj), respectively. 



3.3 Classification of subalgebras 

The classification of subalgebras of is done in the same fashion as for the . The noniden- 
tical adjoint actions involving basis elements of the algebra are exhausted by the following 
list: 

Ad{e'^*)V = V - edt Ad{e'^)dt = e'dt 

Ad(e-'^(f))P = V + eZ{g + tg') Ad{e'^^3)^Q^ = ^ ^z{g') 

Ad(e"^i ) = J2 cos e + J3 sin e Ad{ef-^^)J3 = J3 cos e + Ji sin e 

Ad(e^^i )J3 = -J2sm£ + Js cos e Ad{e'^^)Ji = - J3 sin e + Ji cos e 

Ad{e'^')Ji= cos £ + ^2 sine Ad{e'^' ) Z (g) = Z{g), g = g{t - e) 

Ad{e'^^)J2 = -Ji sine +J2 cose Ad{e''^ ) Z (g) = Z{g), ~g = g{e~H). 

Similar to the case of the /3BVE, we start with the most general form of an infinitesimal 
generator v = a/jP + atdt + a\J\ + 02^/2 + asiTs + Z{g). In the same manner, by acting with 
the adjoint actions given above, we can determine the following list of conjugacy inequivalent 
one- dimensional subalgebras of S^: 

{V + aJi), {dt + aJi), {J^+Z{g)), {Z{g)), (10) 

where a £ M and a G {—1,0, 1} for the first and second cases, respectively. Unlike the case of 
the /9BVE, there is no discrete symmetry allowing placement of additional restrictions on the 
values of a. There are equivalence relations within the two last families of subalgebras, generated 
by adjoint actions of the scaling transformations, time translations, and within the last family, 
changes of algebra bases. 

Using the same procedure as described in the second part of section [2^ we find the following 
list of conjugacy inequivalent two-dimensional subalgebras of (l8|): 

{V + aJi,dt), {V,Ji+Ziat~^)), {V + aJi,Z{\t\^)), 
{dt,Ji+Z{c)), {dt + cJi,Z{e'')), {J,+Zig'),Zig^)), {Z{g'), Z{g^)) , 

where a, 6 G M, c G {—1,0,1}; c G {—1,0,1} if c = 0. There are additional equivalence 
relations within the last two series of subalgebras, generated by adjoint actions of the scale 
transformations, time translations and changes in the algebra bases. 



3.4 Group-invariant reduction with one-dimensional subalgebras 

Based on the above classification of one-dimensional algebras, below we present the correspond- 
ing list of reduced differential equations obtained from the sBVE here. Again, p, q denote the 
new independent variables, while v = v{p, q) is the new dependent variable. 



9 



1 {V + aJi) ijj = t p = X — alnt, q = H 

w + awp - VpWg + VgWp = 0, w := j^Vpp + ((1 - q'^)vg)g 

2 [dt + aJi) ip = V p = X — at, q = /I 

-{aq + v)qWp + {aq + v)pWq = 0, w := j^i^pp + ((1 - q'^)vq)q 

3 {Ji+Z{g)) iP = v + g{t)X p = t, q = 

Wp + gwg = 0, w:= {{I - q'^)vg)q 

4 {Z{g)) No group- invariant reduction is possible in this case 

All Lie symmetries of the reduced equations of Cases 1 and 2 are induced by Lie symmetries of 
the sBVE. This is why further Lie reductions of these cases give no new solutions in comparison 
to Lie reductions with respect to two-dimensional subalgebras. 

We now give some examples for solutions obtained upon using the above ansatze: 

Case 2 includes the well-known Rossby-Haurwitz wave solutions. To show this, we construct 
a class of exact solutions upon using invariance of the sBVE under the algebra {dt + aJi). 
In particular, the corresponding reduced vorticity equation implies that w = F, where is a 
function of v + aq. Hence, we have 

Y^-^Vpp + {1 - q^)vgg - 2qvg = F. (11) 

Eqn. (1111) is, in general, a nonlinear Poisson equation in spherical coordinates. To obtain the 
Rossby-Haurwitz wave solution from this equation, we set F = c{v + aq), c = const, that is, we 
make a homogeneous linear ansatz for F. Separation of the variables gives the ansatz 

vip, q) = Ae'^PPJTiQ) + Be^^PQ^iq) - 

C -r Z 

with A,B = const, where P^{q) and Q^{q) are the associated Legendre functions of the first 
and second kind, respectively, and the degree n is given by 

n = 1(^1^-1). (12) 

For the sake of brevity, we now set B = 0. Reverting to the original variables and employing 
transformation ([9]) to map the solution of the sBVE with vanishing rotation to a solution of the 
sBVE with rotation, we find 

^P{t, X, ^i) = - ^ + ^^M- (13) 

c + I 

To derive pure wave solutions, we require a = Q.{c + 2)/c, which, upon inserting in (|13p and 
considering ()12p . allows us to arrive at the well-known phase relation for a single Rossby- 
Haurwitz wave (e.g. Refs. mill EH): 

20, 

Cphase := a - = -, — TT- (14) 

n[n + 1) 

Since the integer constant m is arbitrary and only linear cases of eqn. (jlip are considered, we 
may extend the solution p3p by superposition of single solutions with m ranging from —n to n. 
This recovers — upon using (jl4p — the classical ansatz for the stream function of Rossby-Haurwitz 
waves. 
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Moreover, note that the class of solutions that may be obtained upon employing symmetry 
methods is again much wider than that obtained upon using the usual ansatz for the stream 
function. In fact, it can be seen that Rossby-Haurwitz waves correspond to particular simple 
solutions of the reduced spherical vorticity equation (jlip. but there is an infinite class of other 
solutions invariant under the same generator dt + aJi . 

Case 3 is completely integrable by quadratures. The general solution for the stream function 
reads 



3.5 Group-invariant reduction with two-dimensional subalgebras 

As was discussed in section 12.51 it is, in general, not necessary to investigate reductions with 
the complete set of two-dimensional inequivalent subalgebras. Namely, if some of the equations 
obtained from one-dimensional reduction are completely integrable, we can avoid the compu- 
tation of reduction with two-dimensional subalgebras if these algebras contain the generators 
that enabled the complete integration in the first place. For the sBVE, case 3 is integrable and 
case 4 does not allow to compute classical group-invariant solutions. Hence, we again have only 
the reduction in one two-dimensional subalgebra, which is not trivial in view of the reductions 
based on one-dimensional subalgebras: (P -|- aJi,dt). Note, however, that it is not possible to 
use this subalgebra in the case a = for a classical Lie reduction since no proper ansatz for 
ip can be constructed. Rather, it can only be used for the construction of partially invariant 
solutions. If a 7^ 0, this subalgebra leads to invariant solutions that are obtainable as particular 
cases of reduction with the algebra {dt). The corresponding ansatz ■0 = e^^v{fi) reduces the 
sBVE to the equation vw^ — v^w = 0, where b = —1/a, w := b^{l — fj?)~^v + ((1 — iJ?)Vfj,)f^. The 
reduced equation implies the following linear constraint between v and w. w = Cv, where C is 
an arbitrary constant, i.e., we have the equation 



which is integrable in terms of Legendre functions. 
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